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EQUATIONS FOR VAPOR PRESSURES AND LATENT 
HEATS, INCLUDING APPROXIMATE EQUATIONS 
FOR SOLID COMPOUNDS CONTAINING 
A GASEOUS COMPONENT.* 


By Louis J. GILLESPIE. 


After early papers by Kirchoff! and Rankine? the subject of vapor 
pressure equations was thoroughly treated by Gibbs,*® subject to 
the provisions that the vapor is an ideal gas of invariable specific 
heat, and the density and the specific heat at constant pressure of 
the liquid or solids are constant. Of the many vapor pressure 
equations that have been published we will mention here only that 
of Nernst,4 which was applied to solids, and which assumes the vapor 
to be an ideal gas, and the specific heats to be functions of the tempera- 
ture (usually in series expansion). Van der Waals and also van Laar 
have discussed the vapor pressure equation’ in the light of the van der 
Waals equation of state. 

An equation derived by Kirchoff and by Gibbs and discussed by 
Gibbs for either solids or liquids is 


log p = A — Blog T — C/T (1) 


where A, B, and C are constants. 

Keyes and Hara® found that log p is linear in 1/7 for the silver- 
silver oxide equilibrium up to 388 atmospheres and pointed out that 
the equation deduced by Pollitzer from the previous measurements of 
Lewis, with the aid of the Nernst heat theorem in connection with 
series expansions of specific heats, failed to represent the data. The 
failure of the equation to serve as an interpolation formula for the 
range covered by the measurements of Lewis is evident: it may be 
that it is not due entirely to the series expansion of specific heats 
at very low temperatures, now known to be erroneous, but possibly 
in part due to failure to apply corrections to the ideal gas law. It 
has often been found that experimental data on decomposition 
pressures may be satisfactorily correlated when B in equation (1) 
is put equal to zero, which corresponds in the derivation of (1) to 
the assumption that the increase of heat capacity in the evaporation 
is zero, and it may be that there is some advantage in either keeping 
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both the term B log 7 and the corrections to the ideal gas law or 
rejecting both together. 

We shall here present the general vapor pressure equation, containing 
no assumptions beyond the two laws of thermodynamics, and then 
develop its application by means of simplifying assumptions, with 
especial reference to the decomposition pressures of solids. 


2. DERIVATION OF THE GENERAL VAPOR PRESSURE EQUATION. 


Let us consider a system in two compartments, separated by a 
membrane which we may at will imagine as solid and impermeable, 
or freely permeable to gas. Let a solid and its compound with a gas, 
for example, BaCl, and BaCl.-8NHs3, be present in one compartment, 
filling it to the exclusion of any gas, and let the pure gas which may 
be evolved by the compound, for example, ammonia, be present alone 
in the other compartment. When the membrane is impermeable we 
may change the condition of the gas by changing independently the 
temperature 7 and the pressure p. The change of chemical potential 
of the gas v. will be given by 


dy = — sdT+vdp (2) 


where s is the entropy and v the volume of a mole of the gas, by Gibbs’ 


equation 97. It will be convenient to write a similar equation for the 
chemical potential of the gaseous component in the compartment 
containing nothing but the solids: 


dy = —3dT+ dp (3) 


where § and @ will refer to the gaseous component in the solids and 
may be called the partial molal entropy and the partial molal volume 
of the gaseous component in the solids. The next four paragraphs 
deal with the exactness of this. | 

The solids, considered collectively in this way as if they were a 
single phase, may be called the condensate. As already pointed out,’ 
“the condensate has then a potential, a partial molal volume, etc., 
with respect to a given component, and the Gibbs theory applies to 
these quantities just as to the usual quantities distributed among the 
individual phases, the arguments not being changed.” This was 
justified by the consideration that there is no logical necessity for 
always dividing up the extensive thermodynamic properties like the 
energy, the thermodynamic potential, the volume, etc., among the 
individual phases, but they may for special purposes be divided up, 
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one portion to the gas, and the rest to the condensed phases considered 
collectively. The reasoning, which was merely outlined, involves 
some difficulties, which will be avoided if we take advantage of Gibbs’ 
equation 121. 

The chemical equation for the formation of the compound may be 
written: 


1/8 BaCl + NH; = 1/8 BaCl-8 NH3 (type reaction) (4) 
vS+G=ySG (general reaction) (5) 


where S represents the chemical symbol for the solid, G that for the 
gas, and SG that for the compound of solid with gas. Writing 
Us, 4, and wsqg for the potentials of solid, gas, and compound, re- 
pectively, Gibbs’ equation 121 becomes 


Vtg = (6) 


Now each potential in equation (6) may be regarded as a function of 
the two variables 7’ and p. We have by differentiation 


(du/dp) 7 = v(dusg/dp) — (dus/dp) (7) 
(du/dp) = v(vgg — Us) (8) 


where vg and vgg are the volume of a mole of solid and compound 
respectively. The last term in equation (8) is the increase in volume 
of the solids attending the reaction represented by equation (5). It 
may be regarded as the partial molal volume of the gaseous component 
in the condensate, since if V be the total volume of the condensate, 
n, the number of moles of gaseous component and nz the number of 
moles of solid S it contains, the partial volume of the gaseous com- 
ponent would be defined as (0V'/dn,), where the variation is at con- 
stant 7, p, and m, and by reason of equation (5) we should have 


(OV = v(%gg — Vg) (9) 


This is therefore, by equation (8), the value of the quantity written 

in equation (3) as 3 for brevity, as the potential y. is the same regard- 

less of the choice of the other components, as shown by Gibbs. 
Similarly it may easily be seen that 


(dy/dT), = — v(8gq — 8s) = — 3 (10) 


It will be observed that the concept of the potential of ammonia 
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in a compound such as BaCl.-8NH; is inconvenient since the partial 
volume of ammonia in it is hardly a fit subject of discussion, as the 
composition of a chemical compound is not susceptible of the continu- 
ous variation indicated by d0E/dn, or 0V/dm. The partial volume 
of ammonia in the condensate is however physically interpretable 
through equation (9). The potential of ammonia in the condensate 
is interpretable through equation (6). The difficulty due to the law 
of definite proportions is similarly removed in the case of other 
partial thermodynamic quantities, such as heat capacity. Such 
considerations apply to the solids of any univariant system, and are 
useful when the system considered has a variance greater than one, 
provided the extra variance is due to the presence of other gases 
which are inert, not reacting with the solids. 

We may now proceed to derive a vapor pressure equation from 
equations (2) and (3) by integration and elimination of the potentials. 
Since du is a perfect differential, the increase of potential Ay of the 
gas in passing from a standard state defined by the values 7) and jx» 
to the variable state given by 7’ and p may be found by integrating 
equation (2) along any path we please, and it is convenient to in- 
tegrate first at the constant pressure pp to the values 7’, po, and then 
at constant temperature 7 to the values 7’, p. The state given by 
7) and py is arbitrary, and will usually be taken at the lowest experi- 
mental temperature and vapor pressure under consideration. - The 
variation of entropy at constant pressure is given by 


(ds/dT), = C,/T (11) 


Hence, if s 9 is the molal entropy of the gas at 7, po, the entropy at 
T, po is 


T 
=f Cid In T + (12) 


Ts 


where C, is the constant pressure specific heat at po and 7. The 
increase of potential Ay for the gas is therefore 


T T 
Ay. = -f Cid \n 7 + 
To To 
p 
+ (v— RT/p) dp + RT In (p/p) (13) 
Do 


Here we have added and subtracted the integral of 27’ dp/p from the 
right-hand side in order to display the term 7 In p, which occurs 


pt 
bay 
: 
mis 
y 
ry 
Fe 
"ig 
rot 
i> 
=) 


EQUATIONS FOR VAPOR PRESSURES AND LATENT HEATS 157 


in any vapor pressure equation as soon as a convenient equation of 
state is chosen. 
For the solids, 


T 
¥ = v(8sg — V(¢sg — dln T + (14) 


by equation (10) and the general equation (11) applied to each 
solid phase, where % is the difference v(ssg — sg) at To, po. Let 
v(e¢sg — ¢s) be represented by ¢, which may be called the molal 
specific heat at constant pressure of the gaseous component in the con- 
densate. Then the increase Ay. for the solids is 


T T 
é dln r+ialart (15) 
To To Po 


Now along the vapor pressure curve y is the same for both gas and 
solid and hence also Ay; furthermore, 7 and p are the same for both, 
including the particular values 7) and po. Elimination of Ay. from 
equations (13) and (15) gives the vapor pressure equation: 


T T 
RT In p = RT In po +f f (C, — @) dln TdT 
To 


p 
+ (89 — — -f (v— RT/p— (16) 
Do 


which may be transformed into 


Jr, (C, — @) din TdT 


Xo 
+ (o— RT/p— (17) 
by virtue of the thermodynamic relation: 
ho = To (80 — 4) (18) 


where Apo is’ the latent heat of vaporization of 1 mole of the gaseous 
component at 79, po. In these equations (C,, — @) is for the constant 
pressure po, and (v — R7'/p — d) is ‘or the temperature 7’, which is 
held constant for the last integration. 

Kquation (17) is the general vapor pressure equation. It is exact, 
depending only on the two laws of thermodynamics. It holds not 
only for solid compounds which evolve a gas, but also for liquids, if 
¢, ¥, and 0 are understood to be the specific heat, the entropy and the 
volume of 1 mole of liquid. 
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3. DISCUSSION OF THE INTEGRALS. 


The value of the single integral of equations (16) and (17) depends 
evidently on both 7 and p. But this integral proves to be of but 
small influence on the calculated value of In p, compared with the 
influence of the term \o/R7’, so that an approximate value of p serves 
to evaluate the integral closely enough to give a very good value of 
In p. It is for this reason that the addition of a quantity equal to 
zero in obtaining equation (13) was significant. 

The single integral of equation (17) may be discussed as regards 
the two parts (v — RT/p) dp and tdp. As regards the second, 3 is 
the volume of the gaseous component in the liquid or in the solids at 
the temperature 7’ and the integration extends from po to p, so that 
the condensed phase would from a physical standpoint be super- 
saturated with vapor all the way up to the upper limit of integration. 
Theoretically this is of no consequence, as a perfect differential (du) 
may be integrated along any path, irrespective of equilibrium con- 
ditions. Practically, the compressibility might be determined if 
necessary at higher pressures and extrapolated into the desired 
range with doubtless sufficient accuracy. 

The integral of (v — RT/p) dp may be calculated from the Beattie- 
Bridgeman equation of state in its original pressure form,® but this 
requires integration by parts, yielding functions of the volume which 
would necessitate solving the equation for the volume by the method 
of successive approximations. 

Beattie® has given another form of the equation, the volume form, 
which is 

RT/p+ B/RT + yp/R?T? + bp?/ RST? (19) 
where 


B = RTBo — Ao — Re/T? 
y = — RTBob + Aga — RBoe/T? (20) 
6 = RBobe/ T? 
which permits an easy integration of (ve — R7/p) dp in a series of 
powers of p: 
(0 — RD /p) dp = — po) 


6 
(p? — yo?) + (py? — (21) 


2R?T? 
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Constants for 14 gases (including air) have been given for the pressure 
equation.!° Constants for the volume equation have not yet been 
_ determined from the pressure data. Beattie has shown however* 
- that constants as determined for the pressure equation may be used 
in the volume equation with good results for gases above their critical 
temperatures. The results were not so good for gases near the 
critical temperature and also near the saturation line, although they 
were considerably better than those obtained with the ideal gas law. 
The double integral of equation (17) is a function of one variable 7’. 
If C, and @ are representable by a series expansion, which will usually 
be the case except at very low temperatures, an expression will 
hold such as 
(22) 


Substitution into the double integral ./2 of equation (17) and integra- 
tion between the limits 75 and 7 gives 


Jo = AT In T — 75 


7,3 
atain Ty + BTy + |r- Ty) 


3 


2-3 


B D ry. ‘ ‘ 
t+ — + — Te") + — Toft) (28) 
4. APPLICATION OF THE VAPOR PRESSURE EQUATION 
To EXPERIMENTAL DATA 


When the specific heats, volumes, and the latent heat are sufficiently 
known, equation (17) furnishes of course a test for the thermodynamic 
consistency of the data. 

When only the specific heats and the volumes are sufficiently 
known, Xo, the latent heat at 79, po, can be found by plotting for 
instance (In p — Jo/RT' + J;/RT), where ./; indicates the single 
integral and J. the double integral of equation (17), against 
(1/7) — 1/7), Xo being found from the slope of the best straight line. 
The constant pp will have been determined from a preliminary smooth- 
ing of the data. 

The case of liquids has not been studied to determine to what 
degree simplification by approximation may safely be carried. In 
this case much is demanded of a vapor pressure equation, on account 
of the high precision attainable in the experimental measurement of 
the vapor pressure of liquids. 
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In the case of solids which evolve a gas, the relations are more 
favorable. The precision attainable in measurements of decomposi- 
tion pressures is not so high as in the case of liquids, probably on 
account of strains in the solids, which are not quite of negligible 
thermodynamic effect. The decomposition pressure of solids at any 
temperature is always less than the normal vapor pressure of the 
liquid, and often much less. Thus there is less demand on the 
equation of state for the gas, and equation (21) will often suffice 
when the constants used are those determined for the pressure form 
of the equation of state. Other approximations will usually be 
permissible, and are made below. 

The quantity 0 is known for many solids for a particular temperature 
and pressure. It does not vary much in a series of similar compounds, 
and when it is unknown for some compound, the value for a similar 
compound may be used. It may be considered a constant with 
respect to temperature and pressure changes. 

As respects the specific heats, the quantity (C, — ¢) will at present 
have to be regarded as a constant, at least in most cases, and as an 
unknown adjustable constant in many cases, because of the lack of 
data for calculating @. 

It appears inadvisable to retain terms expressing variations of 6 or ¢ 
unless these are known from independent data. If they are kept and 
treated as adjustable constants, a better representation of data may 
be expected, but the resulting equation is likely to represent errors 
which ought to be disregarded, and in such a case not only will these 
constants receive fictitious values, but the values found for the latent 
heat will suffer a loss of accuracy. 

Treating (C,, — ) and das constants, and using the volume equation 
for the gas, the following approximate equation is easily obtained 
especially for decomposition pressures of solid compounds. 

M(C,-—7@)f 1 
log p = log po — 1-7) 


Mof1 Mf B 


— My(p? — po?)/2R*T® — M6(p* — pg?)/BRAT4 (24) 


where $8, y, and 6 are calculated from equations (20) for the given 
value of the variable temperature 7', and M is the modulus of loga- 
rithms, 0.43429 . . , common logarithms being denoted by log. 
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If therefore we put 


+ My (p? — + — (25) 


and 
l T, 


To 0 


a plot of y against a should give a straight line whose slope dy/dx is 
— M(C, — ¢)/R and whose intercept on the axis: x = O,is M)d)/R. 
In this way the data can be treated when (C, — ©) is unknown, and 
the constants be determined. The constant po will have been found 
by a preliminary smoothing of the data, plotting log p against 1/7’. 
Points corresponding to values of 7, p very close to 7, po may be 
found far off the line, and no attention need be paid to them in selecting 
the best line. ; 

The constants (C, — ¢) and » having been determined, equation 
24) may be used to calculate the decomposition pressure correspond- 
ing to the temperature 7. The last three terms may be calculated by 
the use of an approximate value of p, as determined for instance by 
use of the preliminary plot or equation used for obtaining the value 
of po, and then corrected if necessary by the use of the improved 
value of p thus found. The convergence is rapid. 

Equation (24) was applied in this way to the data of Gillespie and 
Lurie! on the decomposition pressures of the system 


BaCl. + BaCk-SNHs, 


which extend from 0° C. to about 50° and to about 9 atmospheres. 
No particular difficulty was encountered. In fact, the last two terms 
of equation (24) never exceeded 0.03 per cent of the pressure, a 
negligible matter in the case of an equilibrium involving solids. The 
quantity @ was found to be about 16 calories per mole of ammonia in 
the solids, which may be compared with 11 cal/mole, the value found 
from available data” for AICl; + AICl-6 NH, at about 0° C. A 
value of § was known. 

It seems probable that the terms of equation (24) involving p* 
and p® may often if not usually be omitted, giving, with elimination 
of 8 by equations (20) the relatively simple equation: 
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log p = log pp 2 ( 


Myo v— Bo Ao Cc 
This equation is recommended for solids as the approximation 
usually next in order to equation (1) with B = 0. Equation (1) 
does not appear to have been especially successful and it is but little 
simpler to apply it to data, when B and C are not known in advance, 
than equation (27) itself. 


5. THe Latent Heats oF EVAPORATION. 


The latent heat at 7> in joules or calories will be obtained by 
multiplying Xo, obtained as above in say liter-atmospheres (normal), by 
a factor Q, which may be 101.328 for (absolute) joules, or 24.2122 for 
15°-gram-calories.* 

The latent heat \ at another temperature 7’ is equal to T (s — 3) 
and may be found by applying the exact thermodynamic equations 


(28) and (29): 
T 
So + — (dv/dT),dp (28) 
Do 


= % + edinT - (dv/dT) ,dp (29) 


which give the general oan (30): 
T r P 
i= 7. Ao + rf (C, —e)dinT—-T [d(v —v)/dT]|,dp (30) 
0 To Da 


These integrations are to be performed as in the case of the potentials, 
that with respect to 7 being carried out at the pressure py, and that 
with respect to p being carried out at the temperature T. 

Assuming that (C,, — ¢@) and % are constant and using the equation 
of state (19) we obtain the equation 


To, + M (C, log M log Do 
Ao 3¢ Bob Aga 2Boe 


SBobe 


— Po) (31) 
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to accompany equation (24). In the application to the data of 
Gillespie and Lurie the last two terms of equation (31) were negligible, 
and the approximate formula: 


To + (Cr — ¢) log — 


A 3 
o( + — po) (32) 


is recommended for use in connection with equation (27). In equa- 
tions (30) and (31) C,, ¢, and R will be in the same heat units as \ 
and Xo, and Q is the conversion factor mentioned above. 


6. CONTROLS ON CALCULATIONS. 


The numerical work, the applicability of the simplifying assump- 
tions used and the calculation of the latent heats may be controlled 
as follows. 

Let the vapor pressures be calculated for the experimental tempera- 
tures by the equation derived, and designated by p,. Let the cor- 
responding pressures be calculated by the use of the linear relation (33): 


log p = A— C/T (33) 


and designated by »;. Let the values of the deviation: dev = 
log p. — log »; be plotted against 1/7. They must lie on a smooth 
curve, which may now be drawn. Let the slope of the smooth 
deviation curve at the temperature 7 be d dev/d(1/T). Then 


d log p./d(1/T) = d dev/d(1/T) + d log pf/d(1/T) 
= ddev/d(1/T) —C (84) 
An exact form of the Clapeyron equation is 
d log p/d(i/T) = — MXT/[p(» — d)] (35) 


The last member of equation (34) should be equal to the right-hand 
side of (35). The value of this may be found from the value of \ as 
calculated in pv units and from a value of (v— v) computed from the 
equation of state used. This equality may be tested by drawing a 
line through the point 7’, p of the deviation plot with the slope 
calculated from the heat, namely with slope equal to 


C — MXT/[p(v — 


the line so drawn must be tangent to the deviation curve. 


> a 
aie 
hy, 
is 
Sy 


164 GILLESPIE 


A further control is obtained by plotting the deviation: log p’ — 
log p;, where p’ is the observed pressure, on the same deviation plot. 
If these points fall on a smooth curve of their own, different from the 
other smooth curve, the approximations used for smoothing the 
pressures have been too severe, otherwise not. But if better values of 
d dev/d(1/T) can be thus obtained from the observed points, the 
value of the latent heat in heat units may be calculated from these 
values together with the best possible values of (v — %) from the 
relation: 


where C is the same as in equation (33). 


7. SUMMARY. 


By methods similar to that of Gibbs, general thermodynamic 
equations have been derived for the vapor pressure and the latent 
heats. <A difficulty is discussed and removed, which appears if we 
wish to use a partial quantity such as a partial volume for the case 
of chemical compounds, due to the combination of atoms in definite 
proportions. 

Approximations which may be made when dealing with solid 
compounds are discussed, and approximate equations are given for 
the decomposition pressures of solid compounds and their latent 
heats of decomposition. Methods of applying the equations to 
experimental data are given, together with a procedure for controlling 
the accuracy of the calculations. 
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